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1 Introduction 

Photonic crystals are usually viewed as an optical analog of semiconductors 
that modify the properties of light similar to a microscopic atomic lattice that 
creates a semiconductor band-gap for electrons It is therefore believed 
that by replacing relatively slow electrons with photons as the carriers of 
information, the speed and band-width of advanced communication systems 
will be dramatically increased, thus revolutionizing the telecommunication 
industry. Recent fabrication of photonic crystals with a band gap at optical 
wavelengths from 1.35 (im to 1.95 [im makes this promise very realistic Q. 

To employ the high-technology potential of photonic crystals, it is cru- 
cially important to achieve a dynamical tunability of their band gap || . This 
idea can be realized by changing the light intensity in the so-called nonlinear 
photonic crystals, having a periodic modulation of the nonlinear refractive 
index Exploration of nonlinear properties of photonic band-gap (PBG) 
materials is an important direction of research that opens new applications of 
photonic crystals for all-optical signal processing and switching, allowing an 
effective way to create tunable band-gap structures operating entirely with 
light. 

One of the important physical concepts associated with nonlinearity is the 
energy self-trapping and localization. In the linear physics, the idea of localiza- 
tion is always associated with disorder that breaks translational invariance. 
However, during the recent years it was demonstrated that localization can 
occur in the absence of any disorder and solely due to nonlinearity in the 
form of intrinsic localized modes A rigorous proof of the existence of 
time-periodic, spatially localized solutions describing such nonlinear modes 
has been presented for a broad class of Hamiltonian coupled-oscillator non- 
linear lattices but approximate analytical solutions can also be found in 
many other cases, demonstrating a generality of the concept of nonlinear 
localized modes. 

Nonlinear localized modes can be easily identified in numerical molecular- 
dynamics simulations in many different physical models (see, e.g., Ref. B for 
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a review), but only very recently the first experimental observations of spa- 
tially localized nonlinear modes have been reported in mixcd-valcncc tran- 
sition metal complexes |?]], quasi-one-dimensional antiferromagnetic chains 
M, and arrays of Joscphson junctions ||. Importantly, very similar types 
of spatially localized nonlinear modes have been experimentally observed in 
macroscopic mechanical [|l0] and guided-wave optical jllj systems. 

Recent experimental observations of nonlinear localized modes, as well as 
numerous theoretical results, indicate that nonlinearity-induced localization 
and spatially localized modes can be expected in physical systems of very 
different nature. From the viewpoint of possible practical applications, self- 
localized states in optics seem to be the most promising ones; they can lead 
to different types of nonlinear all-optical switching devices where light ma- 
nipulates and controls light itself by varying the input intensity. As a result, 
the study of nonlinear localized modes in photonic structures is expected to 
bring a variety of realistic applications of intrinsic localized modes. 

One of the promising fields where the concept of nonlinear localized modes 
may find practical applications is the physics of photonic crystals [or pho- 
tonic band gap (PBG) materials] - periodic dielectric structures that produce 
many of the same phenomena for photons as the crystalline atomic potential 
does for electrons [[[]. Three-dimensional (3D) photonic crystals for visible 
light have been successfully fabricated only within the past year or two, and 
presently many research groups are working on creating tunable band-gap 
switches and transistors operating entirely with light. The most recent idea 
is to employ nonlinear properties of band-gap materials, thus creating nonlin- 
ear photonic crystals including those where nonlinear susceptibility is periodic 
as well 

Nonlinear photonic crystals (or photonic crystals with embedded nonlin- 
ear impurities) create an ideal environment for the generation and observa- 
tion of nonlinear localized photonic modes. In particular, the existence of such 
modes for the frequencies in the photonic band gaps has been predicted [[l3f 
for 2D and 3D photonic crystals with Kerr nonlinearity. Nonlinear localized 
modes can be also excited at nonlinear interfaces with quadratic nonlinearity 
| fl4| , or along dielectric waveguide structures possessing a nonlinear Kerr-type 
response |15| . 

In this Chapter, we study self-trapping of light and nonlinear localized 
modes in nonlinear photonic crystals and photonic crystal waveguides. For 
simplicity, we consider the case of a 2D photonic crystal with embedded 
nonlinear rods (impurities) and demonstrate that the effective interaction in 
such a structure is nonlocal, so that the nonlinear effects can be described 
by a nontrivial generalization of the nonlinear lattice models that include the 
long-range coupling and nonlocal nonlinearity. We describe several different 
types of nonlinear guided-wave states in photonic crystal waveguides and 
analyse their properties Jig]. Also, we predict the existence of stable nonlinear 
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localized modes (highly localized modes analogous to gap solitons in the 
continuum limit) in the reduced-symmetry nonlinear photonic crystals |l7j]. 

2 Basic Equations 

Let us consider a 2D photonic crystal created by a periodic lattice of parallel, 
infinitely long dielectric rods in air (see Fig. [I]). We assume that the rods are 
parallel to the X3 axis, so that the system is characterized by the dielectric 
constant e(x) — e(x\,X2)- As is well known 0, the photonic crystals of this 
type can possess a complete band gap for the i?-polarized (with the electric 
field E || 333) light propagating in the (x±, a^-plane. The evolution of such a 
light is governed by the scalar wave equation 

V 2 E(x,t)-^d?{e(x)E}=0, (1) 

where V 2 = 9j + d 2 2 and E is the X3 component of E. Taking the electric 
field in the form E(x,t) = e~ luJt E(x, 1 1 cu) , where E(x,t\ui) is a slowly 
varying envelope, i.e. d 2 E(x, t \ u)) -C u>d t E(x, t\u>), Eq. ([!]) reduces to 



V 2 +£(*)(- 

V c 



21 uj dE 



E(x,t\u) ~ -2ie(x)- — . (2) 

£T Ot 



In the stationary case, i.e. when the r.h.s. of Eq. (g) vanishes, this equation 
describes an eigenvalue problem which can be solved, e.g. by the plane waves 
method Q, in the case of a perfect photonic crystal, for which the dielectric 
constant e(x) = e p (x) is a periodic function defined as 

s p (x + Sij) = e p (x) , (3) 

where i and j are arbitrary integers, and 

s i: j =ia\+ ja 2 (4) 

is a linear combination of the lattice vectors a\ and 

For definitcness, we consider the 2D photonic crystal earlier analyses (in 



the linear limit) in Refs. [g_9 20 1 . That is, we assume that cylindrical rods with 
radius tq = 0.18a and dielectric constant £0 = 11-56 form a square lattice 
with the distance a between two neighboring rods, so that a\ = ax\ and 
ct2 = ax 2 - The frequency band structure for this type of 2D photonic crystal 
is shown in Fig. [I] where, using the notations of the solid-state physics, the 
wave dispersion is mapped onto the Brillouin zone of the so-called reciprocal 
lattice that faces are known as r, M, and X. As follows from Fig. [j], there 
exists a large (38%) band gap that extends from the lower cut-off frequency, 
oj = 0.302 x 2-kc/cl, to the upper band-gap frequency, co = 0.443 x 2nc/a. If 
the frequency of a low- intensity light falls into the band gap, the light cannot 
propagate through the photonic crystal and is reflected. 
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Fig. 1. The band-gap structure of the photonic crystal consisting of a square lattice 
of dielectric rods with ro = 0.18a and Eo = 11.56 (the band gaps are shaded). The 
top center inset shows a cross-sectional view of the 2D photonic crystal depicted 
in the right inset. The bottom center inset shows the corresponding Brillouin zone, 
with the irreducible zone shaded. 



3 Defect Modes: The Green Function Approach 

One of the most intriguing properties of photonic band gap crystals is the 
emergence of exponentially localized modes that may appear within the pho- 
tonic band gaps when a defect is embedded into an otherwise perfect photonic 
crystal. The simplest way to create a defect in a 2D photonic crystal is to 
introduce an additional defect rod with the radius and the dielectric con- 
stant Ed{x). In this case, the dielectric constant e(x) can be presented as a 
sum of periodic and defect-induced terms, i.e. 

e(x) = e p (x) +e d (x) , 

and, therefore, Eq. (|^) takes the form 

where the operator 

t=(^) 2 £d{x) + 2ie{x)%± (6) 

is introduced for convenience. Equation (|^) can also be written in the equiv- 
alent integral form 

E(x,t\u) = J d 2 y G{x 1 y\u J )£E{y 1 t\Lu) , (7) 

where G(x, y \ lo) is the Green function defined, in a standard way, as a solu- 
tion of the equation 



E(x,t\u) = -£E(x,t\w) , (5) 



V 2 + (^)% p (x) 



G(x,y\u) = -5(x-y) . (8) 
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Fig. 2. (a) Electric field structure of a linear localized mode supported by a single 
defect rod with radius = 0.1a and £<j = 11.56 in a square-lattice photonic crystal 
with ro = 0.18a and eo = 11.56. The rod positions are indicated by circles and the 
amplitude of the electric field is indicated by color, (b) Frequency of the defect 
mode as a function of the radius r^: calculated precisely from Eq. (m) {full line with 
circles) and approximately from Eq. (^) [dashed line with triangles). 



General properties of the Green function of a perfect 2D photonic crystal are 
described in more details in Ref. |lq |. Here, we notice that the Green function 
is symmetric, i.e. 

G{x,y\uj) = G(y,x\w) 

and periodic, i.e. 

G(x + s lj ,y + s i0 | uj) = G(x, y | uj) , 

where is defined by Eq. 

The Green function can be calculated by means of the Fourier transform 

/oo 
dt e iut G{x,y,t) (9) 
-OO 

applied to the time-dependent Green function governed by the equation 

[V 2 - e p (x) 9 t 2 ] G(x, y, t) = -S(t) S(x - y) , (10) 

which has been solved by the finite-difference time-domain method |^l) . 

Now that we have calculated the Green function, we can figure out the 
defect states solving Eq. (]?]) directly. For example, Fig. ||(a) shows a defect 
mode created by introducing a single defect rod with the radius r^ — 0.1a 
and dielectric constant — 11.56 into the 2D photonic crystal shown in 
Fig. |l|. Although direct numerical solution of the integral equation (]?]) remains 
possible even in the case of a few defect rods, it becomes severely limited by 
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the current computer facilities as soon as we increase the number of the 
defect rods and start investigation of the line defects (waveguides) and their 
branches. Thus, looking for new approximate numerical techniques which 
could combine reasonable accuracy, flexibility, and power to forecast new 
effects is an issue of the key importance. 



4 Effective discrete equations 

Studying the electric field distribution of the defect mode in Fig. ||(a), one 
can suggest that a reasonably accurate approximation should be provided by 
the assumption that the electric field inside the defect rod remains constant. 
Indeed, let us assume that nonlinear defect rods embedded into a photonic 
crystal are located at the points x m , where m is the index (or a combination 
of two indices in the case of a two-dimensional array of defect rods) introduced 
for explicit numbering of the defect rods. In this case, the correction to the 
dielectric constant is 

e d (x) = {4 0) + \E(x, t\cu)\ 2 }j2 9 ( x - x ™) , (11) 

in 

where 

e(x) = ll> f r \ x \^ rd > (12) 

v ; 1 , for as > r d . K 1 



The second term in Eq. (11) takes into account a contribution due to the 
Kerr nonlinearity (we assume that the electric field is scaled with the nonlin- 
ear susceptibility, x )■ Assuming, as we discussed above, that the electric 
field E(x, t 1 oj) inside the defect rods is almost constant, one can derive, by 
substituting Eq. (|lTJ) into Eq. (0) and averaging over the cross-section of the 
rods (15), an approximate discrete nonlinear equation 

iVg- t E n -E n + J2 Jn-M(e d 0) + \E m \ 2 )E m = , (13) 

m 

for the amplitudes of the electric field E n (t \ oj) = E(x n , 1 1 oj) inside the defect 
rods. The parameter a and the coupling constants 

J„(w) = (^) 2 ^ d 2 y G(x ,x n +y\uj) (14) 

are determined in this case by the Green function G(x,y | oj) of the perfect 
photonic crystal. 

To check the accuracy of the approximation provided by Eq. (|l3|), we 
solved it in the linear limit for the case of a single defect rod. In this case 
Eq. ( |l3| ) is reduced to the equation 

JoM = l/4 0) , (15) 
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Fig. 3. Electric field of the linear guiding mode in a waveguide created by an array 
of the defect rods. The rod positions are indicated by circles and the amplitude of 
the electric field is indicated by color. 



from which one can obtain an estimation for the frequency ojd of the localized 
defect mode. As is seen from Fig. |](b), the mode frequency calculated in the 
framework of this approximation is in a good agreement with that calculated 
directly from Eq. (Q), provided the defect radius is small enough. Even for 
rd = 0.15a an error introduced by the approximation does not exceed 5%. It 
lends a support to the validity of Eq. (|l^) allowing us to use it hereafter for 
studying nonlinear localized modes. 



5 Nonlinear Waveguides in 2D Photonic Crystals 

One of the most promising applications of the PBG structures is a possibil- 
ity to create a novel type of optical waveguides. In conventional waveguides 
such as optical fibers, light is confined by total internal reflection due a dif- 
ference in the refractive indices of the waveguide core and cladding. One of 
the weaknesses of such waveguides is that creating of bends is difficult. Un- 
less the radius of the bend is large compared to the wavelength, much of the 
light will be lost. This is a serious obstacle for creating "integrated optical 
circuits" , since the space required for large-radius bends is unavailable. 

The waveguides based on the PBG materials employ a different physical 
mechanism: the light is guided by a line of coupled defects which possess 
a localized defect mode with frequency inside the band gap. The simplest 
photonic-crystal waveguide can be created by a straight line of defect rods, 
as shown in Fig. ||. Instead of a single localized state of an isolated defect, a 
waveguide supports propagating states (guided modes) with the frequencies 
in a narrow band located inside the band gap of a perfect crystal. Such 
guided modes have a periodical profile along the waveguide, and they decay 
exponentially in the transverse direction, see Fig. ^| That is, photonic crystal 
waveguides operate in a manner similar to resonant cavities, and the light 
with guiding frequencies is forbidden from propagating in the bulk. Because 
of this, when a bend is created in a photonic crystal waveguide, the light 
remains trapped and the only possible problem is that of reflection. However, 
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Fig. 4. The Green function G(xo, xo + y | u>) for the photonic crystal shown in Fig. 
(xq = ai/2 and uj = 0.33 x 2nc/a). 



as was predicted numerically ]19| , f20| and then demonstrated in microwave [ p2[ 
and optical |2^] experiments, it is still possible to get very high transmission 
efficiency for nearly all frequencies inside the gap. 

To employ the high-technology potential of photonic crystal waveguides, 
it is crucially important to achieve a tunability of their transmission prop- 
erties. Nowadays, several approaches have been suggested for this purpose. 
For instance, it has been recently demonstrated both numerically 24 and 
in microwave experiments [^5|, that transmission spectrum of straight and 
sharply bent waveguides in quasiperiodic photonic crystals is rather rich in 
structure and only some frequencies get near perfect transmission. Another 
possibility is creation of the channel drop system on the bases of two parallel 
waveguides coupled by the point defects between them. It has been shown 
[ p6| that high-Q frequency selective complete transfer can occur between such 
waveguides by creating resonant defect states of different symmetry and by 
forcing an accidental degeneracy between them. 

However, being frequency selective, the above mentioned approaches do 
not possess dynamical tunability of the transmission properties. The latter 
idea can be realized by changing the light intensity in the so-called nonlinear 
photonic crystal waveguides pi , created by inserting an additional row of 
rods made from a Kerr-type nonlinear material characterized by the third- 
order nonlinear susceptibility an d the linear dielectric constant e^. For 

defmiteness, we assume that e# = Eq = 11.56 and that the nonlinear defect 
rods are embedded into the photonic crystal along a selected direction Sy , 
so that they are located at the points x m = Xq + m Sjj . As we show below, 
changing the radius r^ of these defect rods and their location Xq in the crystal, 
one can create nonlinear waveguides with quite different properties. 
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Fig. 5. (a) Dispersion relation for the photonic crystal waveguide shown in the 
inset (eo = £d — 11.56, ro = 0.18a, rd = 0.10a). The grey areas are the projected 
band structure of the perfect 2D photonic crystal. The frequencies at the indicated 
points are: uja = 0.378 x 2wc/a and u>b = 0.412 x 2nc/a. (b) Mode power Q(u) of 
the nonlinear mode excited in the corresponding photonic crystal waveguide. The 
right inset gives the dependence Jn(w) calculated at u = 0.37 x 2nc/a. The left 
inset presents the profile of the corresponding nonlinear localized mode. 



As we have already discussed 16 5jJ, the Green function G(x, y \ ui) and, 
consequently, the coupling coefficients J m (w) are usually highly long-ranged 
functions. This can be seen directly from Fig. [| that shows a typical spatial 
profile of the Green function. As a consequence, the coupling coefficients 
J n (<jj) calculated from Eq. ( |l4| ) decrease slowly with the site number n. For 
Soi and s±o directions, the coupling coefficients can be approximated by an 
exponential function as follows 



\M")\ 



Jo(u>) , for n = , 

jjw) e-^W ( for |n| > 1 , 



where the characteristic decay rate a(tu) can be as small as 0.85, depending 
on the values of w, and r^, and it can be even smaller for other types of 
photonic crystals (for instance, for the photonic crystal used in Fig. |lj we 
find J m ~ (— l) m exp(— 0.66m) for m > 2). By this means, Eq. ( fl3| ) is a non- 
trivial long-range generalization of a 2D discrete nonlinear Schrodinger (NLS) 
equation extensively studied during the last decade for different applications 
| j34| . It allows us to draw an analogy between the problem under considera- 
tion and a class of the NLS equations that describe nonlinear excitations in 
quasi-onc-dimensional molecular chains with long-range (e.g. dipole-dipole) 
interaction between the particles and local on-site nonlinearities |27],|28| . For 
such systems, it was shown that the effect of nonlocal interparticle interac- 
tion brings some new features to the properties of nonlinear localized modes 
(in particular, bistability in their spectrum). Moreover, for our model the 
coupling coefficients J n {u>) can be either unstaggered and monotonically de- 
caying, i.e. J n (u>) — \J n (u>)\, or staggered and oscillating from site to site, 




Fig. 6. (a) Dispersion relation for the photonic crystal waveguide shown in the 
inset (eo = £d — 11.56, ro = 0.18a, = 0.10a). The grey areas are the projected 
band structure of the perfect 2D photonic crystal. The frequencies at the indicated 
points are: uja = 0.346 x 1-Kcja and ujb = 0.440 x 2nc/a. (b) Mode power Q(u) of 
the nonlinear mode excited in the corresponding photonic crystal waveguide. Two 
cases are presented: the case of nonlinear rods in a linear photonic crystal (solid 
line) and the case of a completely nonlinear photonic crystal (dashed line). The right 
inset shows the behavior of the coupling coefficients J n (u>) for n > 1 (Jo = 0.045) 
at u) = 0.33 x 2-Kc/a. The left inset shows the profile of the nonlinear mode. 



i.e. J n (u>) = (— 1)"| J n (uj)\. We therefore expect that effective nonlocality in 
both linear and nonlinear terms of Eq. ([l3]) may also bring similar new fea- 
tures into the properties of nonlinear localized modes excited in the photonic 
crystal waveguides. 



5.1 Staggered and unstaggered localized modes 

As can be seen from the structure of the Green function presented in Fig. |], 
the case of monotonically varying coefficients J n (uj) can occur for the waveg- 
uide oriented in the Soi direction with Xq = ai/2. In this case, the frequency 
of a linear guided mode, that can be excited in such a waveguide, takes a 
minimum value at k — [see Fig. p(a)], and the corresponding nonlinear 
mode is expected to be unstaggered. 

We have solved Eq. fll5| ) numerically and found that nonlinearity can lead 
to the existence of guided modes localized in both directions, i.e. in the direc- 
tion perpendicular to the waveguide, due to the guiding properties of a chan- 
nel waveguide created by defect rods, and in the direction of the waveguide, 
due to the nonlinearity- induced self-trapping effect. Such nonlinear modes 
exist with the frequencies below the frequency of the linear guided mode of 
the waveguide, i.e. below the frequency uja hi Fig. ^|(a), and are indeed un- 
staggered, with the bell-shaped profile along the waveguide direction shown 
in the left inset of Fig. mb). 



Nonlinear Localized Modes in Photonic Crystals 



11 




Fig. 7. (a) Dispersion relation for the photonic crystal waveguide shown in the 
inset (eo = £d — 11.56, ro = 0.18a, rd = 0.10a). The grey areas are the projected 
band structure of the perfect 2D photonic crystal. The frequencies at the indicated 
points are: u)a = 0.352 x 2nc/a, lub = 0.371 x 2nc/a, and u>c = 0.376 x 2vrc/a (at 
k — 0.217 x 2n/a). (b) Mode power Q(lj) of the nonlinear mode excited in the 
corresponding photonic crystal waveguide. The right inset shows the behavior of 
the coupling coefficients J n (u>) for n > 1 (Jo = 0.068) at uj — 0.345 x 2irc/a. The 
left inset shows the profile of the corresponding nonlinear mode. 



The 2D nonlinear modes localized in both dimensions can be characterized 
by the mode power which we define, by analogy with the NLS equation, as 

Q = Y J \En\ 2 - (16) 

n 

This power is closely related to the energy of the electric field in the 2D 
photonic crystal accumulated in the nonlinear mode. In Fig. ^](b) we plot 
the dependence of Q on frequency, for the waveguide geometry shown in Fig. 
|(a). 

As can be seen from the Green function shown in Fig. [|, the case of stag- 
gered coupling coefficients J n (to) can be obtained for the waveguide oriented 
in the s w direction with x = a\/2. In this case, the frequency dependence 
of the linear guided mode of the waveguide takes the minimum at k = ir/a 
[see Fig. §(a)]. Accordingly, the nonlinear guided mode localized along the 
direction of the waveguide is expected to exist with the frequency below the 
lowest frequency loa of the linear guided mode, with a staggered profile. The 
longitudinal profile of such a 2D nonlinear localized mode is shown in the left 
inset in Fig. [|(b), together with the dependence of the mode power Q on the 
frequency (solid curve), which in this case is again monotonic. 

The results presented above are obtained for linear photonic crystals with 
nonlinear waveguides created by a row of defect rods. However, we have 
carried out the same analysis for the general case of a nonlinear photonic 
crystal that is created by rods of different size but made of the same nonlinear 
material. Importantly, we have found relatively small difference in all the 
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Fig. 8. Examples of the (a) symmetric and (b) antisymmetric localized modes. 
The rod positions are indicated by circles and the amplitude of the electric field is 
indicated by color [red, for positive values, and blue, for negative values]; (c) Power 
Q vs. frequency dependencies calculated for two modes of different symmetry in the 
photonic crystal waveguide shown in Fig. ^ 



results presented above provided nonlinearity is weak. In particular, for the 
photonic crystal waveguide shown in Fig. ^(a), the results for linear and 
nonlinear photonic crystals are very close. Indeed, for the mode power Q the 
results corresponding to a nonlinear photonic crystal are shown in Fig. ^(b) 
by a dashed curve, and for Q < 20 this curve almost coincides with the solid 
curve corresponding to the case of a nonlinear waveguide embedded into a 
2D linear photonic crystal. 



5.2 Instability of nonlinear localized modes 

Let us now consider the waveguide created by a row of defect rods which 
are located at the points Xo — (ai + 0,2) /2, along a straight line in cither 
the S10 or S01 directions. The results for this case are presented in Figs. [7] |[ 
The coupling coefficients J n are described by a slowly decaying staggered 
function of the site number n, so that the effective interaction decays on the 
scale larger than in the two cases considered above. 
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It is remarkable that, similar to the NLS models with long-range dispersive 
interactions p7[p^ |, we find a non-monotonic behavior of the mode power 
Q(u>) for this type of nonlinear photonic crystal waveguides: specifically, Q(u>) 
increases in the frequency interval 0.344 < (wa/27rc) < 0.347 [shaded in 
Fig. ||(c)]. One can expect that, similar to the results earlier obtained for the 
nonlocal NLS models |^7],^8j, the nonlinear localized modes in this interval 
are unstable and will eventually decay or transform into the modes of higher 
or lower frequency j?9). What counts is that there is an interval of mode 
power in which two stable nonlinear localized modes of different widths do 
coexist. Since the mode power is closely related to the mode energy, one can 
expect that the mode energy is also non-monotonic function of lu. Such a 
phenomenon is known as bistability, and in the problem under consideration 
it occurs as a direct manifestation of the nonlocality of the effective (linear 
and nonlinear) interaction between the defect rod sites. 

Being interested in the mobility of the nonlinear localized modes we have 
investigated, in addition to the symmetric modes shown in the left inset in 
Fig. R(b) and in Fig. ^(a), also the antisymmetric localized modes shown in 
Fig. g(b). Our calculations show that the power Q(ui) of the antisymmetric 
modes always (for all values of uj and all types of waveguides) exceeds that for 
symmetric ones [see, e.g., Fig. ||(c)]. Thus, antisymmetric modes are expected 
to be unstable and they should transform into a lower-energy symmetric 
modes. 

In fact, the difference between the power of antisymmetric and symmetric 
modes determines the Peierls-Nabarro barrier which should be overtaken for 
realizing the mobility of a nonlinear localized mode. One can see in Fig. 
||(c) that the Peierls-Nabarro barrier is negligible for 0.347 < (yiajlixc) < 
0.352 and thus such localized modes should be mobile. However, the Peierls- 
Nabarro barrier becomes sufficiently large for highly localized modes with 
lo < 0.344 x 2irc/a and, as a consequence, such modes should be immobile. 
Hence, the bistability phenomenon in the photonic crystal waveguides of the 
type depicted in Figs. opens up fresh opportunities p8| for switching 
between immobile localized modes (used for the energy storage) and mobile 
localized modes (used for the energy transport). 

The foregoing discussions on the mode mobility, based on the qualitative 
picture of the Peierls-Nabarro barrier, have been established for the discrete 
one- dimensional arrays. It is clear that the two-dimensional geometry of pho- 
tonic crystals under consideration will bring new features into this picture. 
However, all these issues are still open and would require a further analysis. 

6 Self- Trapping of Light in a Reduced-Symmetry 2D 
Nonlinear Photonic Crystal 

A low-intensity light cannot propagate through a photonic crystal if the light 
frequency falls into a band gap. However, it has been recently suggested Q 
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Reduced-symmetry 
square lattice 




Fig. 9. Band-gap structure of the reduced-symmetry photonic crystal with ro = 
0.1a, Td = 0.05a, and e = 11.4 for both types of rods. Full lines are calculated 
by the MIT Photonic-Bands program j3l) whereas dashed line is found from the 
effective discrete model. The top center inset shows a cross-sectional view of the 
2D photonic crystal depicted in the right inset. The bottom center inset shows the 
corresponding Brillouin zone. 



that in the case of a 2D periodic medium with a Kerr-type nonlinear material, 
high-intensity light with frequency inside the gap can propagate in the form of 
finite energy solitary waves - 2D gap solitons. These solitary waves were found 
to be stable JT^| , but the conclusion was based on the coupled-mode equations 
valid for a weak modulation of the dielectric constant e(x). However, in real 
photonic crystals the modulation of e(x) is comparable to its average value. 
Thus, the results of Ref. jl3| have a limited applicability to the properties of 
localized modes in realistic photonic crystals. 

More specifically, the coupled-mode equations are valid if and only if the 
band gap A is vanishingly small, i. e. A ~ A 2 where A is an effective am- 
plitude of the mode, that is a small parameter in the multi-scale asymptotic 
expansions J32| . If we apply this model to describe nonlinear modes in a wider 
gap (see, e.g., discussions in Ref. J3^|), we obtain a 2D nonlinear Schrodinger 
(NLS) equation known to possess no stable localized solutions. Moreover, the 
2D localized modes described by the coupled-mode equations are expected 
to possess an oscillatory instability recently discovered for a broad class of 
coupled-mode Thirring-like equations p3| ). Thus, it is clear that, if nonlinear 
localized modes do exist in realistic PBG materials, their stability should be 
associated with different physical mechanisms not accounted for by simplified 
continuum models. 

In this Section we follow Ref. JD] and study the properties of nonlinear 
localized modes in a 2D photonic crystal composed of two types of circular 
rods: the rods of radius ro made from a linear dielectric material and placed 
at the corners of a square lattice with the lattice spacing a, and the rods 
of radius r^ made from a nonlinear dielectric material and placed at the 
center of each unit cell (see right inset in Fig. ||). Recently, such photonic 
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Fig. 10. Coupling coefficients J n ,m(w) for the photonic crystal depicted in Fig. 
^ (the contribution of the coefficient Jo,o = 0.039 is not shown). The frequency 
uj = 0.4456 falls into the first band gap. 



crystals of reduced symmetry have attracted considerable interest because of 
their ability to possess larger absolute band gaps [^p| . The band-gap structure 
of the reduced-symmetry photonic crystal is shown in Fig. As is seen, it 
possesses two band gaps, first of which extends from u> — 0.426 x 2nc/a to 
lo = 0.453 x 2nc/a. 

The reduced-symmetry "diatomic" photonic crystal shown in Fig. ^| can 
be considered as a square lattice of the "nonlinear defect rods" of small radius 
Td ("I'd < r o) embedded into the ordinary single-rod photonic crystal formed 
by a square lattice of rods of larger radius ro in air. The positions of the 
defect rods can then be described by the vectors x n ^ m — na\ +ma,2, where 
a\ and a-i are the primitive lattice vectors of the 2D photonic crystal. Here, in 
contrast to the photonic crystal waveguides discussed in the previous section, 
the nonlinear defect rods are characterized by two integer indices, n and m. 
However, it is straightforward to extend Eq. ([l3]) and write an approximate 
2D discrete nonlinear equation 

ia—E^ m -E n ^ + Y,Jn-k^ m -i(to){e { d 0) + \E k d\ 2 )E k ^=0, (17) 
k,l 

for the amplitudes of the electric field E n:Tn (t \ lo) = E(x n , m , 1 1 lo) inside the 
defect rods. We have checked the accuracy of the approximation provided by 



Eq. (17) solving it in the linear limit, in order to find the band-gap struc- 
ture associated with linear stationary mode. Since the coupling coefficients 
Jn,m{^) in the photonic crystal depicted in Fig. ^| are highly long- ranged 
functions (see Fig. |l^), one should take into account the interaction between 
at least 10 neighbors to reach accurate results. As is seen from Fig. [|, in this 
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Fig. 11. Top (left) and 3D (right) views of a nonlinear localized mode in the first 
band gap of 2D photonic crystal depicted in Fig. 0. 



case the frequencies of the linear modes (depicted by a dashed line, with a 
minimum at u> — 0.446 x 2irc/a) calculated from Eq. ( p^ ) are in a good agree- 
ment with those calculated directly from Eq. (Q). It lends a support to the 
validity of Eq. (|l7]) and allows us to use it for studying nonlinear properties. 

Stationary nonlinear modes described by Eq. ( |l7| ) are found numerically 
by the Newton-Raphson iteration scheme. We reveal the existence of a con- 
tinuous family of such modes, and a typical example [smoothed by continuous 
optimization for Eq. (|^)] of nonlinear localized mode is shown in Fig. [H]. In 
Fig. |l^, we plot the dependence of the mode power 

QH = ^K, m | 2 , (18) 

n,m 

on the frequency u) for the photonic crystal shown in Fig. |[ As we have 
already discussed, this dependence represents a very important characteristic 
of nonlinear localized modes which allows to determine their stability by 
means of the Vakhitov-Kolokolov stability criterion: dQ/duj > for unstable 
modes (this criterion has been extended |35| to 2D NLS models). 

As is well known |35|]36[| , in the 2D discrete cubic NLS equation, only 
high-amplitude localized modes are stable, whereas no stable modes exist in 
the continuum limit. For our model, the high- amplitude modes are also stable 
(see inset in Fig. |l2|), but they are not accessible under realistic conditions: 
To excite such modes one should increase the refractive index at the mode 
center in more than 2 times. Thus, for realistic conditions and relatively small 
values of > only low-amplitude localized modes become a subject of much 
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Fig. 12. Power Q vs. frequency uj for the 2D nonlinear localized modes in the 
photonic crystal of Fig. [] with two different e^ ' . Solid lines - stable modes, dashed 
lines - unstable modes. Insets show typical profiles of stable modes, and an enlarged 
part of the power dependence. Grey areas show the lower and upper bands of 
delocalized modes surrounding the band gap. 



interest since they can be excited in experiment. However, such modes in 
unbounded 2D NLS models are always unstable and either collapse or spread 
out Q. In fact, they can be stabilized by some external forces (e.g., due to 
interactions with boundaries or disorder |37J ) , but in this case the excitations 
are pinned and cannot be used for energy or signal transfer. 

Here we reveal that, in a sharp contrast to the 2D discrete NLS models 
discussed earlier in various applications, the low-amplitude localized modes 
of Eq. ( jl7|) can be stabilized due to nonlinear long-range dispersion inherent 
to the photonic crystals. It should be emphasized that such stabilization does 
not occur in the models with only linear long-range dispersion P4| . In order 
to gain a better insight into the stabilization mechanism, we have carried out 
the studies of Eq. (|l3| ) for the exponentially decaying coupling coefficients 
Jn,m- Our results show that the most important factor which determines 
stability of the low-amplitude localized modes is a ratio of the coefficients 
at the local nonlinearity (~ Jo,o) an d the nonlinear dispersion (~ Jo.i)- If 
the coupling coefficients J rhm decrease with the distances n and m rapidly, 

the low-amplitude modes of Eq. ( p~7| ) with eji ' = 11.4 are essentially stable 
for Jq,o/Jq,i < 13. However, this estimation is usually lowered because the 
stabilization is favored by the presence of long-range interactions. 

It should be mentioned that the stabilization of low-amplitude 2D local- 
ized modes is not inherent to all types of nonlinear photonic crystals. On the 
contrary, the photonic crystals must be carefully designed to support stable 
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low- amplitude nonlinear modes. For example, in the photonic crystal consid- 
ered above such modes are stable at least for 11 < < 12, however they 
become unstable for > 12 (see Fig. |l2|). The stability of these modes can 
also be controlled by varying r^, ro, or eo- 

7 Concluding Remarks 

Exploration of nonlinear properties of PBG materials may open new im- 
portant application of photonic crystals for all-optical signal processing and 
switching, allowing an effective way to create tunable band-gap structures op- 
crating entirely with light. Nonlinear photonic crystals, and nonlinear waveg- 
uides created in the photonic structures with a periodically modulated dielec- 
tric constant, create an ideal environment for the generation and observation 
of nonlinear localized modes. 

As follows from our results, nonlinear localized modes can be excited in 
photonic crystal waveguides of different geometry. For several geometries of 
2D waveguides, we have demonstrated that such modes are described by 
a new type of nonlinear lattice models that include long-range interaction 
and effectively nonlocal nonlinear response. It is expected that the general 
features of nonlinear guided modes described here will be preserved in other 
types of photonic crystal waveguides. Additionally, similar types of nonlinear 
localized modes are expected in photonic crystal fibers |58| consisting of a 
periodic air-hole lattice that runs along the length of the fiber, provided the 
fiber core is made of a highly nonlinear material (see, e.g., Ref. [p9|). 

Experimental observation of nonlinear photonic localized modes would 
require not only the use of photonic materials with a relatively large nonlinear 
refractive index (such as AlGaAs waveguide PBG structures or polymer 
PBG crystals MJ], but also a control of the group- velocity dispersion and 
band-gap parameters. The latter can be achieved by employing the surface 
coupling technique that is able to provide coupling to specific points of the 
dispersion curve, opening up a very straightforward way to access nonlinear 
effects. 
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